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THERMAL STRUCTURE OF THE MARTIAN 
ATMOSPHERE DURING THE DISSIPATION 
OF THE DUST STORM OF 1971 

Abstract 

The secular variation of the thermal structure of the Martian atmosphere 
during the dissipation phase of the 1971 dust storm is examined, using tem- 
peratures obtained by the infrared spectroscopy investigation on Mariner 9. 

For the latitude i-ange -20° to -30°, the mean temperature at the 2mb level is 
found to decrease from approximately 220K in mid-December 1971 to about 
190K by June 1972 while for the 0.3mb level a decrease from 203K to 160K is 
observed. Over the same period, the amplitude of the diurnal temperature wave 
also decreased. Assuming a simplified radiative heating model, the dust optical 
depth is found to decrease approximately exponentially with an e-folding time of 
about 60 days at both the 0.3 and 2mb levels, Stokes -Cunningham settling alone 
cannot account for this behavior. Sedimentation models which include both 
gravitational settling and vertical mixing are developed in an effort to explain 
the time evolution of the dust. Within the framework of a model which assumes 
an effective vertical diffusivity K independent of height, a mean dust particle 
diameter of ~ 2 pm is inferred. To provide the necessary vertical mixing, 

K > 10 7 cm 2 sec " 1 is required in the low^er atmosphere. 
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THERMAL STRUCTURE OF THE MARTIAN 
ATMOSPHERE DURING THE DISSIPATION 
OF THE DUST STORM OF 1971 


Introduction 

The Mariner 9 spacecraft achieved orbit about Mars 14 Nov. 1971 during a 
planet -wide dust storm of an extraordinary magnitude. In television pictures 
returned by the spacecraft during the first few weeks of the mission, the surface 
of the planet was found to be almost totally obscured. During the later part of 
December, a clearing trend began with surface features becoming increasingly 
well defined. By mid-January, surface features were visible almost everywhere, 
even at large viewing angles (Leovy, et al., 1972). 

Data obtained from the infrared spectroscopy investigation on Mariner 9 
have been used to obtain atmospheric temperature profiles during the dust storm 
and its subsequent decay (Hanel, et al., 1972a; 1972b; Conrath, et al., 1973). 

Data acquired during the storm revealed that the atmosphere was in a highly 
perturbed thermal state compared to that expected for a clear atmosphere on 
the basis of theoretical models (Gierasch and Goody, 1969). This suggested that 
absorption of solar energy by the dust was strongly affecting the atmospheric 
energy balance, (Gierasch and Goody, 1972). A diurnal temperature wave with 
a peak-to-peak amplitude of ~30K was observed at all atmospheric levels up to 
at least 40km, and mean temperatures were substantially higher than those 
predicted for a clear atmosphere. These observations have been used to deduce 
the properties of the strong thermally driven atmospheric tides which existed 
during the storm (Leovy, et al., 1973; Pirraglia and Conrath, 1974). 
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In the present paper, the secular behavior of the atmospheric temperature 
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d following the storm is examined; The data are 


interpreted in terms of a simplified atmospheric heating model from which a 
time dependent dust optical depth is deduced. In an effort to explain the behavior 
of the optical depth as a function of time, a model describing the dynamics of 
the dust settling process is invoked. This model contains both gravitational 
settling and vertical mixing of the dust particles and is used to estimate particle 
size and effective vertical diffusivity. 


Data Analysis 

Temperature profiles were calculated using measurements from selected 
spectral intervals within the 6^7 cm -1 CC> 2 absorption band. The techniques 
employed are similar to those developed for use in the terrestrial meteorological 
satellite program (see for example Conrath and Revah, 1972). Temperatures 
were obtained from the surface to ~0.1mb with a vertical resolution on the order 
of half of a pressure scale height. 

The temperature data were organized according to latitude and local time. 

In order to define the secular trend, the data prior to orbit 220 (the first 110 days 
of the mission) were averaged over intervals of approximately six days. Examples 
of these averages are shown in Figures 1 and 2 for the 0.3mb and 2mb atmos- 
pheric levels respectively in the latitude interval -20° to -30°. Following orbit 
220, a much smaller quantity of data was acquired, necessitating the use of 
individual data points rather than averages. Nevertheless, the continuing trend 
of the data appears to be well defined. 
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Examination of Figures 1 and 2 reveals that the diurnal mean temperature 
at both levels decreased monotonically approaching ~160K at 0.3mb and ~ 90 K 
at 2mb. Comparison of the behavior of data taken in early morning with that in 
late afternoon permits an estimate of the amplitude of the diurnal wave which is 
also found to decrease with time. This suggests that the presence of dust is 
responsible for enhanced atmospheric heating. To relate the temporal behavior 
of the atmospheric thermal structure directly to the behavior of the dust, it is 
necessary to first develop an atmospheric heating model. 


Radiative Heating Model 

A highly simplified radiative heating model will now be considered. It is 
assumed that the heat equation can be written 


3T _ (T-T e ) Q d 
t r + C p p 


( 1 ) 


where T is atmospheric temperature which is a function of height z and time t , 
T is the radiative-convective equilibrium temperature which the atmosphere 

e 

would assume in the absence of dust, p is atmospheric density, C p is specific 
heat at constant pressure, t r is the radiative-convective time constant of the 
* atmosphere and Q D is atmospheric heating due to the dust. It is assumed that 
the perturbed thermal state of the atmosphere is produced by direct deposition 
of solar energy into the atmosphere due to the presence of the dust. The 
volumetric heating is approximated by 


Qd 


dF 

dz 


(2) 
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where F is the solar flux reaching level z which can be expressed in the form 

f }jl S e‘ r/fl ; 0 < n < 1 

F=< (3) 

Lo ; -1 < (jl £ 0 . 


S is the solar flux at the top of the atmosphere, r is the extinction optical depth 
and jj. is the cosine of the solar zenith angle. It is assumed that a fraction a of 
the solar energy removed from the beam goes into heating the atmosphere and 
solar energy reflected from the surface has been neglected. A similar heating 
model has been used previously by Leovy, et al. (1972). 

With the above assumptions , the heating rate per gram of atmosphere re- 
quired in (1) can be written 


Op 

^ C p 


with the definition 


G(z, t) e’ T < z » t > / ^< t > ; 0 < yLi< 1 

0 ; -1 < /u. < 0 


(4) 


G (z, t ) = f ^g( t ) q (z, t) (5) 

P 

where q{z, t) is the specific concentration of the dust (number of dust particles 
per gram of atmosphere) and <x is the extinction cross section for the dust. Use 
has been made of the relationship between optical depth and specific concentration 
which can be written 


CO 

T(Z, t) = cr q p d z , 


( 6 ) 
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The time dependent parameters of the model vary on two distinct time 
scales; The cosine of the solar zenith angle at colatitude 9 and time t is 
given by 

^(t) = b(t) cos nt + c(t) (7) 


where 0 is the planetary rotation rate, b(t) = sin 6> s (t) sin 9 and c(t) = cos 
Q (t) cos e with 6 s (t) being the (seasonally varying) latitude of the subsolar 
point. Thus, b, c, s, q, and r are all slowly varying on a time scale of one day. 
This suggests that the secular variations can be separated from the diurnal 
fluctuations by assuming 


and 


T(t) -T e = A 0 (t) + 


L 


[A n (t) cos nflt + B n (t) sin nflt] 


( 8 ) 


exp [— t(z, t) />(t)] = a 0 (t) + L 3 n cos n t ( ) 

n 

where the coefficients in the expansions are slowly varying on a time scale of 
D” 1 . Then 


h>(t) = - 


77 


r. i d 4- 

l b(t) cos + c (t ) J 




r < 10 ) 


(t ) f exp [~ r ^. z -’— — 1 cos n <p 

" J Q L ‘-'V*-/ ^ y' w v " / 


dcp 


J 


where 4> =fit and 4>. = cos" 1 (-c/b). The coefficients employed in the subsequent 
calculations were evaluated from (10) using numerical integration. 
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Substituting (8) into (1) and using (4), (5) and (9), one obtains 


aA o A o 


dt 

+ t 

r 

dA 

A 

n 

n 

dt 

+ t f 

dB 

B 

n 

^ n 


B n =G(z,t) a n (t), > 


J 


( 11 ) 


The assumption that the coefficients change slowly on a time scale of fl” 1 , 
along with the assumption that t r is of the order of one day or less allows the 
first terms on the left hand sides of each of equations (11) to be neglected. The 
coefficients can then be obtained and the solution to (1) written in the form 

T(z,t) = T(z, t ) + 2 ^ AT n( t ) COS ( nQt - (12 > 

n 

where 


T(z, t) = T e (z) + t r G(z,t) a 0 (t) 


(13) 


AT n (z, t ) = t f G(z, t ) a n (t )//l + n 2 Q 2 t 2 


The equilibrium temperature T e for each level is taken to be the apparent 
asymptotic value approached by the temperature data in Figures 1 and 2. The 
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remaining parameters of the model are obtained by fitting the mean temperature 
T (t) and the diurnal amplitude AT, (t) to the observed data. A good fit is found 
by assuming 

r(z, t) = t( z, 0) e 7 c (14) 

The resulting T(t) (solid line) and T(t) ±ATj(t) (dashed lines) are shown in 
Figures 1 and 2 for t c =60 days for both levels, and t r = 5 x 10 4 sec. It is 
noteworthy that a! good fit is obtained at both levels with the same characteristic 
time t c . 

Dust Dissipation 

Having estimated the temporal behavior of the extinction optical depth for 
solar radiation, it is of interest to consider the implications which can be drawn 
concerning the dust dissipation process. If it is assumed that during the dis- 
sipation phase no new dust particles are injected into the atmosphere and all 
particles reaching the ground remain there, then the dissipation is controlled 
by gravitational settling and vertical mixing. Since the gravitational settling 
rate is dependent on dust particle size, and the vertical mixing is generally 
parameterized in terms of an effective diffusivity, it should be possible in 
principle to obtain information on these parameters. To do this, it is necessary 
to define a model describing the dust dissipation. 

It is assumed that the continuity equation governing the behavior of the 
suspended dust can be written 


3q + 1_ 33 

3t p 3z 


( 15 ) 
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The first term in (16) represents that portion of the total flux due to gravitational 
settling, and the second term is the contribution due to vertical mixing. The 
terminal velocity of the dust particles is denoted by w*. All processes con- 
tributing to the vertical mixing of the dust are parameterized in terms of an 
effective vertical diffusivity K. 

The terminal velocity required in (16) is assumed to be given by the Stokes- 
Cunningham relation for spherical particles (Cunningham, 1910; Green and Lane, 
1957; Ryan, 1964) 



where p p is dust particle density, g is the Martian gravitational acceleration, 
d is particle diameter, 77 is molecular viscosity, and k is molecular mean free 
path. The pressure dependence of k is approximated by 

\(p) = \(p r )^l (18) 

where p r is an arbitrary reference pressure. In subsequent calculations, the 
following numerical values are used: p p = 3 gm cm -3 , g = 370cm sec " 2 , 

V = 1.5 X 10 " 4 gm sec -1 cm' 1 , and k (25mb) = 2.2 x 10' 4 cm. For the range 
of pressure and particle sizes considered here 2k/d » 1 so to good approxi- 
mation w* can be taken as proportional to p _1 d. 
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It is convenient to define a height-dependent characteristic settling time 


t. (O = 


H 

w*(z ) 


(19) 


where H is the pressure scale height. In the small particle limit where w* cc p" 1 , 
t s (z) is the time required for the particle to fall from infinity to level z. In the 
large particle limit, w* is independent of z, and t s is the time required for the 
particle to fall through one scale height. Values of t s for the 0.3, 2, and 5mb 
levels are plotted as functions of particle diameter in Figure 3. In this calcu- 
lation and in all subsequent calculations, it is assumed that H = 10km. 

Consider first the conditions which existed during the storm prior to the 
onset of dissipation. In this case, it seems reasonable to assume a steady state 
(i.e., 3q/B t = 0), and therefore 3 = constant. If it is assumed that dust was 
injected into the atmosphere as rapidly as it was removed at the surface, then 
the net flux at the surface must vanish and 3 = 0 everywhere giving 


w*q_K^ = 0. (20) 

o z 

Before proceeding further, it is necessary to make an assumption concerning 
the height dependence of K. At this point it will be assumed that K is independent 
of height, and the effects of this assumption on the final results will be examined 
later. With the particle terminal velocity approximated by 

w* = - w* e z/H (21) 
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where -w* is the value of w* at z = 0, the solution of (20) is 


q(z) = q 0 ex P [i - exp (^-')]} (22) 

where q 0 is the value of q at z = 0, and 

Hw; 

v= — • (23) 

If a characteristic diffusion time is defined at t D = H 2 /K f then v = t D /t s (0) 
where from (19) t s (0) is the characteristic settling time at z = 0. 

Examples of steady-state dust distributions calculated from (22) for selected 
values of v are given in Figure 4. As the mixing is increased, and hence v 
become smaller, and the dust extends to greater heights. In the limit of infinitely 
rapid mixing, v -» 0 and the specific concentration of the dust becomes in- 
dependent of height. 

It is convenient to define an effective ’'top" of the dust layer Z t at the point 
of inflection in the distribution (22) which can be written Z t = -H 1 n v. This is 
the level at which the settling time t s (z) becomes equal to the diffusion time tjj. 
To maintain a dust distribution up to a given level Z t , the value of K required 
is dependent on the particle diameter. This dependence, calculated using (23) 
along with the definition of z t , is shown in Figure 5 for several values of Z t . 
Television pictures of the limb of the planet during the storm indicate dust up 
to at least 50km (Leovy, et al., 1972). This suggests Z t = 50km corresponding 
to v = .007 in (22). 
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Turning now to the problem of modeling the dissipative phase of the storm, 
time dependent solutions of (15) will be sought. It will be assumed that at time 
t = 0, the mechanism for injecting dust into the atmosphere is switched off, and 
all particles reaching the surface remain there. The distribution for v = .007, 
which is shown as the broken curve in Figure 4, will be taken as the initial value 
of q at t = o. 

First consider the case when there is no mixing (K = 0), and (15) becomes 


Bq 1 B 
Bt + p Bz 


O' * p q ) = o . 


(24) 


Using (21) and letting £ = p/p Q = e" z/H , the solution to (24) can be written 



(25) 


where f(£) is an arbitrary function to be determined from the initial conditions. 
Using (22) and requiring that q be non-negative, 


r 


qo exp{,[l-exp (- f . t 1 /t (o) )]} ‘ 


q(<f. t) = < 


(26) 


1 ° 


; £< 


t 


The optical depth, calculated from (6) , is 
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r 

CT q 0 Po e ' 


t 

J- ^ . 

H 

r ^ i 

t 

- >- » 

L t s (0) _ 

^2 

— i 

/'“N 

o 

'w' 

on 

< 

t 

t s (0) 




(27) 


0 


; £< 


t 

M°) 


where E 2 (x) is the exponential integral of second order defined as 


E 2 (x) = 




dy. 


Figure 6 shows the optical depths at the 0.3 and 2mb levels as functions of time 
calculated from (27) with v = .007. The optical depth decreases much more 
rapidly at the upper level than at the lower level. Obviously this does not agree 
with the behavior deduced from the temperature data. Therefore, it must be 
concluded that vertical mixing played a significant role in the behavior of the 
dust during the dissipative phase of the storm as well as during the active phase. 

It is now necessary to obtain solutions to (15) retaining both the gravitational 
settling and diffusive terms in the particle flux. Using (21) and p = p 0 e“ z/H , 

(15) can be written 


Bq _ B 2 q 

Bt' 


+ (ve?'-l) 


(28) 


where the dimensionless indpendent variables are t' = t/t D and z* = z/H. Again 
the initial value is assumed to be given by (22) with v = .007. Note, however, 
that the value of v appearing in (28) can in general have different value than that 
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used in the initial condition since the value of K may have been different during 
the dissipative phase. In addition to the initial value, boundary conditions must 
also be specified. For z -* 00 it is assumed that q o. Specification of the 
lower boundary condition has been discussed in connection with sedimentation 
problems in the terrestrial atmosphere by several authors (Johnstone, et al., 
1949; Monin, 1959; Calder, 1961). Here it will be assumed that the diffusive 
portion of the flux must vanish at z = 0. Since a constant value of K has been 
assumed, this implies that (Bq/3z) z=0 = 0. Thus, dust particles are ultimately 
removed from an atmospheric column by pure gravitational settling out of the 
bottom of the column onto the planetary surface where they remain. 

Solutions to (28) were obtained numerically for the initial value and boundary 
values discussed above. An implicit method of solution described by Richtmyer 
and Morton (1967) was employed. In practice q was required to vanish at 
z = 100km. However, the exact choice of the level for which q vanishes is not 
critical provided it is higher than the effective "top" of the initial distribution. 

Calculations were made for a range of values of v . An example of the time 
evolution of q with v - .01 is shown in Figure 7. Optical depths at the 0.3 and 
2mb levels obtained from q using (6) are show'n in Figure 8 as functions of time. 
Examination of these curves indicates that the optical depth at the upper level 
decreases significantly more rapidly than that at the lower level for all values 
of v greater than about .01. Thus, in order to obtain agreement with the optical 
depths inferred from the atmospheric temperature data it is necessary to choose 
v < .01. 

From Figure 8d, the time dependence of the optical depths is essentially 
exponential with an e-folding time of approximately t/t s (o) = 1.2. For t = 60 
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days as inferred from the temperature data, t s (o) = 50 days. Assuming an 
atmospheric pressure at the surface of 5mb, this settling time corresponds to 
a dust particle diameter of 2.3 (Figure 3). From a knowledge of the settling 
time and the limit established on v, it is found from (23) that K £ 2 x lo 7 cm 2 -' 
sec -1 . With a knowledge of the particle diameter, it is also possible to estimate 
the value of K required during the active phase of the storm to maintain dust up 
to ~50km. From Figure 5 with d ~ 2/xm, it is found that K ~ 10 7 cm 2 sec -1 
is needed. Thus, the model considered here requires similar values of K during 
both the active and dissipative phases of the storm. 

The assumption of constant K has resulted in large values for this parameter, 
corresponding to strong mixing at all atmospheric levels up to at least 50km. To 
investigate the dependence of the results on the height variation of K, calculations 
were made for two additional models. Results obtained with K(z) = K 0 exp(z/2H) 
are shown in Figure 9 for w*H/K 0 = .05; larger values of this parameter result 
in a more rapid decay in r at the upper level. From the figure, the e-folding 
time is t/t s (o) y 0.9 which for t = 60 days results in d = 1.8 /xm and K 0 = 4 x 10 6 
cm 2 sec -1 . This particular form of K(z) is of interest since it is representative 
of mixing induced by vertically propagating gravity waves (Cf. Wofsy and McElroy, 
1973). As an extreme case of height dependence, solutions were also obtained 
with K(z) = K Q e l/H with the results shown in Figure 10 for w* H/K 0 = 0.2. 
Again requiring an e-folding time of 60 days, it is found that K 0 = 9 x 10 s cm 2 
see -1 and d y 1.4 ^m. 

The diffusivity profiles obtained with the three different models are 
summarized in Figure 11. Also shown is the profile required for the diffusion 
time to be equal to the gravitational settling time at all levels. It is not 



possible to derive a unique height dependence for K on the basis of the data 
analyzed here. However, the basic picture emerging from the calculations can 
be summarized as follows. The secular dependence of the optical depths as 
inferred from atmospheric temperature measurements can be explained in 
terms of models in which the dust is strongly mixed up to ~ 50 km. The 
residence time of the dust is determined primarily by the gravitational settling 
rate at the bottom of an atmospheric column. For this reason, the particle 
size inferred from the residence time does not depend strongly on the height 
profile of K. To maintain a strongly mixed state within an atmospheric column, 
the characteristic diffusion time must be of the order of or less than the gravi- 
tational settling time at each level. The values of K shown represent the mini- 
mum amount of mixing required for each model and therefore represent lower 
limits. Between the surface and 50 km, the mean values of K for the two height 
dependent models are of the same order as the value of K obtained with the 
height independent model. 

The particle size estimate of d ~ 2 obtained here is comparable with 
most of the previous estimates which have been made using other techniques. 
Moroz and Ksanfomaliti (1972), using infrared photometry at 0.7 and 1.4 am 
from Mars 3, estimate the particle size to be about lp.m while Conrath, et al. 
(1973) find a mean particle radius between 1 and 10 from the characteristics 
of the mineral dust features in the infrared spectra. Pang and Hord (1973) 
estimate the mean radius to be about 2 /xm based on ultra violet spectrometer 
data from Mariner 9. Leovy, et al. (1972) infer somewhat larger particles with 
mean radii in excess of 10/xm, based on television photometry from Mariner 9. 
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Recently, Hartmann and Price (1974) have examined the time dependence of the 
contrast of surface features observed in the Mariner 9 television pictures, 
obtaining particle sizes £ 5 pm in the upper atmosphere and £ 10pm in the 
lowest layers. 

Conclusions 

Atmospheric temperatures derived from the infrared spectroscopy investi- 
gation on Mariner 9 during the dissipation of the dust storm of 1971 reveal the 
manner in which the atmosphere relaxed from the perturbed thermal state w’hich 
existed during the storm. It is found that the observed secular changes in mean 
temperatures and diurnal amplitudes can be well fit with a highly simplified 
radiative heating model. Within the framework of this model, it is possible to 
estimate the time dependence of the optical depth of the dust above the 0.3mb 
and 2mb levels. An approximately exponential behavior with the same charac- 
teristic time (60 days) is found for each level. 

Since Stokes-Cunningham settling alone would result in clearing at a more 
rapid rate above the 0„3mb level than above the 2mb level, it is concluded that 
vertical mixing must have played a significant role during the dissipation phase. 
Applications of dust dissipation models which include both gravitational settling 
and vertical mixing, indicate that the mean diameter of the particles responsible 
for the heating was approximately 2 pm. It is not possible to uniquely establish 
the height dependence of the vertical diffusivity K from the data considered here. 
With the assumption of a height independent diffusivity, it is found that K £ 10 7 
cm 2 sec “ 1 is required to maintain sufficient vertical mixing to produce the 
observed behavior of the optical depths during the dissipation of the storm. 
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A similar value of K would have been required during the active phase of the 
storm to maintain a steady state distribution of 2 /an size particles up to 50km. 

ACKNOWLEDGEMENTS 

The author is indebted to J. A. Pirraglia and P. J. Gierasch for useful 
discussions of the subject matter of this paper. The Martian atmospheric 
temperatures were obtained in collaboration with V. G. Kunde, J. C. Pearl and 
R. A. Hanel, Principle Investigator of the Mariner 9 Infrared Spectroscopy 
Investigation, who also provided a critical review of the manuscript. H. W. Sagges 
provided the computer programming necessary for carrying out the numerical 
calculations. 

REFERENCES 

C alder, K. L. (1961). Atmospheric diffusion of particulate material considered 
as a boundary value problem. Journal of Meteorology 18 , 413-416. 

Conrath, B. J., and Revah, I. (1972). A review of nonstatistical techniques for 
the estimation of vertical atmospheric structure from remote infrared 
measurements. In "Proceedings of a Workshop on the Mathematics of 
Profile Inversion" (L. Colin, Ed.), NASA TMX-62150. 

Conrath, B., Curran, R., Hanel, R., Kunde, V., Maguire, W., Pearl, J., Pirraglia, 
J., and Welker, J. (1973). Atmospheric and surface properties of Mars 
obtained by infrared spectroscopy on Mariner 9. J. Geophys. Res. 78, 
4267-4278. 

Cunningham, E. (1910). On the velocity of steady fall of spherical particles 
through fluid medium. Proc. Roy. Soc. London (Series A) 83, 357-365. 


17 




Gierasch, P. and Goody, R. (1968). A study of the thermal and dynamical 

structure of the Martian lower atmosphere. Planet. Space Sci. 16, 615-646. 

Gierasch, P., and Goody, R. (1972). The effect of the dust on the temperature 
of the Martian atmosphere. J. Atmos. Sci. 29 , 400-402. 

Green, H. L., and Lane, W. R. (1957). ’’Particulate Clouds”, E. & F. N. Spon 
Ltd., London. 

Hanel, R. A., Conrath, B. J., Hovis, W. A., Kunde, V. G., Low T man, P. D., 

Pearl, J. C., Prabhakara, C., and Schalchman, B. (1972a). Infrared 
spectroscopy experiment on the Mariner 9 mission: Preliminary results. 
Science 175 , 305-308. 

Hanel, R., Conrath, B., Hovis, W., Kunde, V., Lowman, P., Maguire, W., 

Pearl, J., Pirraglia, J., Prabhakara, C., and Schlachman, B. (1972b), 
Investigation of the Martian environment by infrared spectroscopy on 
Mariner 9. Icarus 17 , 423-442. 

Hartmann, W. K. , and Price, M. J. (1974). Mars: Clearing of the 1971 dust 
storm. Icarus 21 , 28-34. 

Johnstone, H, F,, Winsche, W. E., and Smith, L, W., (1949). The dispersion and 
deposition of aerosols. Chem. Rev. 44 , 353-371. 

Leovy, C. B., Briggs, G. A., Young, A. T., Smith, B. A., Pollack, J. B., Shipley, 
E. N., and Widey, R. L. (1972). The Martian atmosphere: Mariner 9 
television experiment progress report. Icarus 17, 373-393. 

Leovy, C. B., Zureck, R. W., and Pollack, J. B. (1973). Mechanisms for Mars 
dust storms. J. Atmos. Sci. 30, 749-762. 


18 



Monin, A. S. (1959). On the boundary condition on the earth surface for 
diffusing pollution, in "Advances in Geophysics" (Vol. 6), pp. 435-436. 
Academic Press, New York. 

Moroz, V. I., and Ksanfomaliti, L. V. (1972). Preliminary results of astro- 
physical observations of Mars from Mars 3. Icarus 17, 408-422, 

Pang, K., and Hord, C. W. (1973). Mariner 9 Ultraviolet spectrometer experi- 
ment: 1971 Mars’ dust storm. Icarus 18, 481-488. 

Pirraglia, J. A., and Conrath, B. J. (1974). Martian tidal pressure and wind 
fields obtained from the Mariner 9 infrared spectroscopy experiment. 

J. Atmos, Sci. 31 , 318-329. 

Richtmyer , R. D., and Morton, K. W. (1967). "Difference Methods for Initial- 
Value Problems" (2nd edition), pp. 198-201. Interscience, New York. 

Ryan, J. A. (1964). Notes on the Martian yellow clouds. J. Geophys. Res. 69, 
3759-3770. 

Wofsy, S. C., and McElroy, M. B. (1973). On vertical mixing in the upper 
stratosphere and lower mesosphere. J. Geophys. Res. 78 , 2619-2624. 


19 



oon 


220 
210 
200 
£ 190 

s 

K 180 
170 


UJ 

or 


< 

ce 


160 

150 


* *\ 


* 6 TO 7 HRS * 
o 7 TO 8 HRS * 
a 8 TO 9 HRS « 

* 9 TO 10 HRS « 

* 10 TO 11 HRS « 


13 TO 14 HRS 

14 TO 15 HRS 

15 TO 16 HRS 

16 TO 17 HRS 

17 TO 18 HRS 



X 


100 


x 


x 


x 


200 


500 


600 


300 400 

ORBIT NUMBER 

NOV | DEC | JAN | FEB | MAR | APR | MAY | JUN | JUL [ AUG | SEPTj OCT 


EARTH DATE 


Figure 1. Secular variation of atmospheric temperatures at the 0.3 mb level for 
the latitude interval -20° to -30°. The abscissa is Mariner 9 orbit number with 
the corresponding earth date also shown. Data points prior to orbit 210 repre- 
sent averages over approximately 6 days while data points after orbit 210 repre- 
sent individual measurements. The data have been orgainized by Marti an local 
time as indicated in the key to the various symbols. The solid line represents 
the mean temperature and the broken lines the limit of the diurnal excursion 
calculated using a heating model described in the ,text. 
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EARTH DATE 

Figure 2. Secular variation of atmospheric temperatures at the 2 mb level for 
the latitude interval -20° to -30°. The abscissa is Mariner 9 orbit number with 
the corresponding earth date also shown. Data points prior to orbit 210 repre- 
sent averages over approximately 6 days while data points after orbit 210 repre- 
sent individual measurements. The data have been organized by Martian local 
time as indicated in the key to the various symbols. The solid line represents 
the mean temperature and the broken lines the limit of the diurnal excursion 
calculated using a heating model described in the text. 
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PARTICLE DIAMETER (jum) 

Figure 3. Gravitational settling times for dust particles versus particle diam- 
eter for three atmospheric pressure levels. The settling time is defined by 
equation (19) in the text. 
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Figure 4. Vertical dust profiles obtained assuming steady state 
conditions and a constant vertical diffusivity K. Curves are shown 
for selected values of the parameter v, the ratio of the character- 
istic time to the gravitational settling time at the planetary surface. 
The specific dust concentration q is normalized to its value at the 
surface. 
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Figure 6. Calculated time dependence of optical depths at the 0.3 and 2 mb leve 
assuming gravitational settling only. The optical depth at each level is norma 
.zed to its value at time t = 0. The time scale is normalized to t s (0), the grm 
tational settling time at the surface. 
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Figure 7. Example of the time evolution of the vertical dust profile obtained 
with v - .01 in the dissipation model discussed in the text. The specific dust 
concentration q is normalized to its value at the surface at time t = 0, and the 
times indicated for each curve are expressed in terms of the gravitational 
settling time at the surface. 
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v = 0-1 

Figure 8. Calculated time dependence of optical depths at the 0.3 and 2 mb levels assuming 
both gravitational settling and vertical mixing. The optical depth at each level is normalized 
to the gravitational settling at the surface. Results are shown for selected values of v , the 
ratio of the characteristic diffusion time to the gravitational settling time at the surface: 
(a) v = 0.1, (b) v = 0.05; (c) v = 0.02; (d) v = 0.01. A height independent vertical diffusivity 

was assumed. 





Figure 8. (continued). 
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(d) v= 0.01 

Figure 8. (concluded). 
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Figure 10. Calculated time dependence of optical depths at the 0.3 and 2 mb 
levels assuming a vertical diffusivity K = Ko exp{z/H). In this example the 
ratio of the characteristic diffusion time to the gravitational settling time at 
the surface is 0.2. 
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Figure 11. Comparison of the three models of K employed to obtain the results 
shown in Figures 8, 9, and 10. The broken line indicates the value of K required 
for the characteristic diffusion time to be equal to the gravitational settling time 
at each level for 2 /x particles. 
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